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ABSTRACT 

We propose a scenario for the formation of the Main Belt in which asteroids incor- 
porated icy particles formed in the outer Solar Nebula. We calculate the composition 
of icy planetesimals formed beyond a heliocentric distance of 5 AU in the nebula by 
assuming that the abundances of all elements, in particular that of oxygen, are solar. 
As a result, we show that ices formed in the outer Solar Nebula are composed of a mix 
of clathrate hydrates, hydrates formed above 50 K and pure condensates produced at 
lower temperatures. We then consider the inward migration of solids initially produced 
in the outer Solar Nebula and show that a significant fraction may have drifted to the 
current position of the Main Belt without encountering temperature and pressure con- 
ditions high enough to vaporize the ices they contain. We propose that, through the 
detection and identification of initially buried ices revealed by recent impacts on the 
surfaces of asteroids, it could be possible to infer the thermodynamic conditions that 
were present within the Solar Nebula during the accretion of these bodies, and during 
the inward migration of icy planetesimals. We also investigate the potential influence 
that the incorporation of ices in asteroids may have on their porosities and densities. 
In particular, we show how the presence of ices reduces the value of the bulk density 
of a given body, and consequently modifies its macro-porosity from that which would 
be expected from a given taxonomic type. 

Key words: solar system: formation - minor planets, asteroids. 



1 INTRODUCTION 

In recent years, some objects within the Main Belt of as- 
teroids have been found to display cometary characteristics 
(Hsieh & Jewitt 2006). Objects such as 133P/Elst-Pizarro, 
P/2005 Ul and 118401 (1999 RE 70 ) occupy orbits that are 
entirely decoupled from Jupiter within the Main Belt, and 
are probably bodies that have undergone a recent collision, 
revealing previously buried volatile material, and leading to 
the observed dusty outgassing. In addition, present-day sur- 
face water ice and possible water sublimation have been re- 
ported on Ceres (Lebofsky et al. 1981; A'Hearn & Feldman 
1992; Vernazza et al. 2005). This is consistent with recent 
Hubble Space Telescope (HST) observations which suggest 
that Ceres' shape is the result of the dwarf planet consisting 
of a rocky core surrounded by an ice-rich mantle (Thomas 



et al. 2005) - an idea in agreement with several thermal evo- 
lution scenarios (McCord & Sotin 2005) that suggest that 
the ice content of the asteroid is between 17% and 27%, 
by mass. These observations are supported by the evidence 
of hydrated minerals in meteorites which provide samples of 
rock from asteroids in the Main Belt. Most of these minerals 
formed as a result of water ice accreting with the chondritic 
meteorite parent bodies, melting, and driving aqueous al- 
teration reactions (Clayton & Mayeda 1996; Jewitt et al. 
2007). 

It seems likely, then, that some objects in the aster- 
oid belt have incorporated significant amounts of water ice 
(and possibly other volatiles) during their formation in the 
early stages of the solar System. These bodies would have 
incorporated icy particle^] coming from the outer nebula 
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1 By icy particles is meant planetesimals composed of a mix of 
ices and rocks. 
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that survived their inward drift due to gas-drag through the 
disk (Mousis & Alibert 2005 - hereafter MA05). The volatile 
fraction incorporated in this manner could vary depending 
on the inward flux of icy planetesimals from the external 
region and the heliocentric location of the asteroid, together 
with the density of the proto-solar nebula. Given that the 
current asteroid belt lies closer to the Sun than the "snow- 
line" , postulated to lie at around 5 AU in the Solar Nebula, 
these results are a little unexpected. In this context, under- 
standing how volatiles were incorporated into the asteroids 
is therefore important, not only for the study of the asteroids 
themselves, but also for our understanding of the processes 
by which the solar system came into being. 

To this end, MA05 studied the possibility of determin- 
ing the nature and composition of the ices which were in- 
corporated into Ceres. They used a time dependant model 
of the Solar Nebula and showed that icy particles of sizes 
between 0.1 and 10 metres could drift from heliocentric dis- 
tances greater than 5 AU to the present location of Ceres 
without encountering temperatures or pressures high enough 
to vaporise the ices within. The authors then suggested that 
ices produced in the outer Solar Nebula were transported in- 
wards to become incorporated in the solids which accreted 
to form Ceres. 

The present work aims to improve upon the calculation 
detailed in MA05, along with expanding the results to in- 
volve the entire asteroid belt, rather than just its largest 
member. In particular, MA05 postulated that all volatiles, 
except CO2EL were trapped by water in the form of hy- 
drates or clathrate hydrates in the outer solar nebula. This 
assumption was supported by the work of Hersant et al. 
(2001) who estimated that Jupiter was formed at tempera- 
tures higher than ~40 - 50 K. The accretion of ices in the 
form of hydrates and clathrate hydrates was thus required 
during the formation of the planet in order to explain the 
volatile enrichments observed in its atmospher^fl (Gautier 
et al. 2001a,b). Indeed, since these ices are usually formed 
at temperatures higher than that reached by the nebula at 
the time of Jupiter's completion, as defined by Hersant et 
al. (2001), they can be incorporated in the planetesimals ac- 
creted by the giant planet during its growth. However, the 
amount of water that would be required in the nebula to trap 
all these volatiles as hydrates and clathrate hydrates exceeds 
that derived from the solar oxygen abundance. Therefore, 
MA05 made the ad hoc hypothesis that oxygen was "overso- 
lar" in the gas-phase in order to provide enough available wa- 
ter in the Solar NebulsQ. Additionally, Hersant et al. (2001) 
only used an evolutionary Solar Nebula model to derive the 



2 CO2 is the only major volatile species which does not form a 
clathrate hydrate in the Solar Nebula because it condenses as a 
pure ice prior to being trapped by water. 

3 The abundances of volatile species in Jupiter's atmosphere have 
been measured using the mass spectrometer on board the Galileo 
probe. These measurements reveal that the giant planet's atmo- 
sphere is enriched by a factor of ~3 in Ar, Kr, Xe, C, N, and S 
compared to the solar abundances (Owen et al. 1999). 

4 The oxygen abundance required for to allow the trapping of 
all volatile species in the form of hydrates or clathrate hydrates is 
~1.9 times the solar abundance, with C02:CO:CH4 = 1:1:1 and 
N2:NH3 = 1:1 (the nominal nebula gas phase ratios used in this 
work) . 



disk's temperature at the time when the mass of Jupiter's 
feeding zone was equal to that of the gas in its current en- 
velope. They thus neglected many important effects such as 
the influence of protoplanet formation on the structure of 
the disk (e.g. Fig. 2 of Alibert et al. 2004). However, recent 
giant planet core-accretion formation models that include 
migration, disk evolution, such as that proposed by Alibert 
et al. (2004), have shown that the disk's temperature can 
be as low as ~ 10 - 20 K at the end of Jovian formation. 
This implies that Jupiter itself can accrete ices during its 
formation that were produced at temperatures lower than 
those required for clathration. As a result, no extra water 
is required in the nebula to allow all the volatile species to 
be trapped in clathrate hydrates, and the oversolar oxygen 
abundance condition in the nebula can be relaxed. 

In Section 2, we calculate the composition of ices pro- 
duced in the outer Solar Nebula under the assumption that 
the abundances of all elements, in particular that of oxy- 
gen, are solar. In Section 3, we consider the inward migra- 
tion of particles produced at various locations in the nebula, 
and at different times. This allows us to examine whether 
some planetesimals formed in the outer Solar Nebula may 
have drifted to the current position of the Main Belt with- 
out encountering temperature and pressure conditions high 
enough to vaporize the ices they contain. In Section 4, we 
examine the uncertainties in the determination of the physi- 
cal properties of asteroids. We also investigate the potential 
influence that the incorporation of ices in these objects may 
have on their porosities and densities. Section 5 is devoted 
to summary and discussion. 



2 COMPOSITION OF ICY PLANETESIMALS 
FORMED IN THE OUTER SOLAR NEBULA 

2.1 Initial gas-phase conditions in the nebula 

In this Section, we aim to estimate the composition of icy 
particles produced in the region of giant planet formation, 
prior to the dissipation of the nebula, which then migrated 
into the forming Main Belt. The calculations have been car- 
ried out in a manner consistent with the formation of Jupiter 
with a realistic primordial volatile composition. This implies 
that the icy planetesimals that drifted inwards to the for- 
mation zone of the asteroids shared the same composition 
as those accreted by proto- Jupiter. 

We then assume that the gas-phase abundances of el- 
ements are solar (e.g. Table [TJ and that O, C, and N exist 
only in the form of H 2 0, CO2, CO, CH 4 , N 2 , and NH 3 . S is 
only present as H2S and other refractory sulfur components 
(Pasek et al. 2005). The gas-phase molecular ratios in the 
Solar Nebula are presumed to derive directly from that in 
interstellar ices (Mousis et al. 2002) and, in some cases, from 
the consideration of catalytic effects that might affect these 
ratios in the nebula. Thus, although the CO:CH4 ratio is 
typically around 10:1 in the ISM (Allamandola et al. 1999), 
this can be revised down in the Solar Nebula gas-phase since 
some additional CH4 may be introduced through conver- 
sions of both CO and CO2, as a result of the presence of 
catalytically active regions in the disk (Kress & Tielens 2001; 
Sekine et al. 2005). We therefore adopt CO:CH 4 = 1:1 in the 
vapor phase of the nebula. Moreover, CO2 should initially 
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Table 1. Gas phase abundances (molar mixing ratio with respect 
to H2) of major species in the Solar Nebula from Lodders (2003) 
with C0 2 :CO:CH 4 = 1:1:1 and N 2 :NH 3 = 1:1. 



Species X 


X/H 2 


species X 


X/H 2 





1.16.10" 3 


N 2 


5.33. 10~ 5 


c 


5.82.10" 4 


NH 3 


5.33. 10~ 5 


N 


1.60.10" 4 


S 


3.66. 10~ 5 


H 2 


5.78.10" 4 


Ar 


8.43. 10" 6 


CO2 


1.94.10" 4 


Kr 


4.54.10" 9 


CO 


1.94.10" 4 


Xe 


4.44.10" 10 


CH 4 


1.94.10" 4 







be present in the gas-phase of the nebula, with C0 2 :CO = 
1:1 - 4:1, a range of values that covers the ISM measure- 
ments (Gibb et al. 2004). The value of the N 2 :NH 3 ratio 
is quite uncertain, although current chemical models of the 
ISM predict that molecular nitrogen should be much more 
abundant than ammonia (Irvine & Knacke 1989). On the 
other hand, the N 2 :NH3 ratio may have been much lower in 
the Solar Nebula since the conversion of N 2 into NH3 can be 
accelerated by the catalytic effect of local Fe grains (Lewis 
& Prinn 1980; Fegley 2000). In all our following calculations, 
we consider N 2 :NH3 = 1:1 as the nominal ratio for these two 
molecules in the nebula gas-phase. 

2.2 Formation sequence of the different ices in the 
outer nebula 

As a result of the adoption of a solar gas-phase oxygen abun- 
dance, we show here that ices formed at distances greater 
than 5 AU in the cooling Solar Nebula are a mix of clathrate 
hydrates and hydrates formed above 50 K, and pure conden- 
sates primarily produced at lower temperatures (but still 
greater than ~20 K). We note that this hypothesis is sup- 
ported by the recent work of Mousis & Marboeuf (2006), 
who showed that the abundances of volatiles observed in 
the envelopes of Jupiter and Saturn could be reproduced by 
using a solar abundance for all elements in the Solar Nebula, 
and that the calculated amount of heavy elements remains in 
agreement with internal structure models. The clathration 
and hydratation processes result from the presence of avail- 
able crystalline water ice at the time of volatile-trapping in 
the Solar Nebula. This statement seems justified when one 
considers that current scenarios of the formation of the Solar 
Nebula suggest that most ices falling from the presolar cloud 
onto the disk were vaporised upon entering the early nebula. 
Following Chick and Cassen (1997), H 2 ice is initially va- 
porized at locations in the inner 30 AU of the Solar Nebula. 
With time, the decrease of temperature and pressure led to 
conditions which allowed the water to re-condense, forming 
microscopic crystalline ices (Kouchi et al. 1994, Mousis et 
al. 2000). Once formed, the different ices agglomerated and 
were incorporated into the growing and drifting particles. 

The process by which volatiles are trapped, illustrated 
in Fig.[T] is calculated using the stability curves of hydrates, 
clathrate hydrates and pure condensates, and the tracks de- 
tailing the evolution of temperature and pressure at helio- 
centric distances of 5, 15 and 20 AU within the solar nebula. 



Table 2. Parameters of the stability curves of the considered 
clathrate hydrates (reproduced from Hersant et al. 2004). Their 
equations are of the form In P = A/T + B, where P and T are the 
partial pressure (bars) and the temperature (K) of the considered 
species, respectively. A is in K and B is dimensionless. 



Species 


A 


B 


CH 4 


-2161.81 


11.1249 


CO 


-1685.54 


10.9946 


N 2 


-1677.62 


11.1919 


NH 3 


-2878.23 


8.00205 


H 2 S 


-3111.02 


11.3801 


Ar 


-1481.78 


9.95523 


Kr 


-1987.5 


9.99046 


Xe 


-2899.18 


11.0354 



These evolutionary tracks are derived from the Q-turbulent 
accretion disk model described in next section. The stability 
curves of hydrates and clathrate hydrates derive from Lunine 
& Stevenson (1985)'s compilation of laboratory data, from 
which data is available at relatively low temperature and 
pressure. Their equations, given by Hersant et al. (2004), 
are of the form In P = A/T + B, where P and T are the 
partial pressure (bars) and the temperature (K) of the con- 
sidered species, respectively. Table[2]shows the values of con- 
stants A and B determined by Hersant et al. (2004) from 
their fits to laboratory measurements. On the other hand, 
the stability curves of pure condensates used in our calcula- 
tions derive from the compilation of laboratory data given 
in the CRC Handbook of chemistry and physics (Lide 2002). 
Their equations are of the form log P — A/T + B, where 
P and T have the same units as mentioned above. Table [3] 
also gives the constants A and B as determined by our own 
fits to laboratory measurements. The cooling curve inter- 
cepts the stability curves of the different ices at some given 
temperature, pressure and surface density conditions (see 
Table 3]). For each ice considered, the domain of stability is 
the region located below its corresponding stability curve. 
The clathration process stops when no more crystalline wa- 
ter ice is available to trap the volatile species. For example, 
if one assumes C0 2 :CO:CH 4 = 1:1:1 and N 2 :NH 3 = 1:1 in 
the gas-phase, then the NH 3 and H 2 S are entirely trapped 
by the water, as hydrates of NH 3 and clathrate hydrates of 
H 2 S, while only approximately half of the CH 4 is trapped 
as a clathrate hydrate in the Solar Nebula. The remaining 
CH 4 , as well as Xe, Kr, CO, Ar and N 2 , whose clathration 
processes normally occur at higher temperatures, remain in 
the vapor phase until the Solar Nebula cools enough to al- 
low the formation of pure condensates (roughly between 20 
and 30 K in the gas-phase conditions assumed for the So- 
lar Nebula; see Fig. [TJ). Note that, during the cooling of the 
Solar Nebula, C0 2 is the only species that crystallizes as a 
pure condensate prior to being trapped by water to form a 
clathrate hydrate. Hence, we assume here that solid C0 2 is 
the only existing condensed form of C0 2 in these environ- 
ments. 



4 Mousis et al. 



Table 3. Parameters of the stability curves of the considered 
pure condensates. Their equations are of the form log P = A/T 
+ B, where P and T are the partial pressure (bars) and the 
temperature (K) of the considered species, respectively. A is in K 
and B is dimensionless. 



Table 4. Temperature T(K), total gas pressure P(bars) and sur- 
face density S(g/cm 2 ) conditions at which different ices form in 
the Solar Nebula at heliocentric distance of 5 AU. The gas-phase 
conditions considered here are C02:CO:CH4 = 1:1:1 and N2:NH3 
= 1:1. 



Species 


A 


B 


Ices 


T(K) 


P(bar) 


S(g/cm 2 ) 


CH 4 


-475.61 


4.2831 


H a O 


157.5 


3.4 X 10~ 7 


677.4 


CO 


-411.24 


5.2426 


NH 3 -H 2 


89.1 


1.9 x 10" 7 


498.0 


co 2 


-1365.9 


7.0248 


H 2 S-5.75H 2 


84.0 


1.8 x 10~ 7 


482.6 


N 2 


-360.07 


4.7459 


co 2 


78.0 


1.6 x 10~ 7 


463.8 


NH 3 


-1565.0 


6.7883 


CH 4 -5.75H 2 


60.9 


1.3 x 10~ 7 


404.3 


H 2 S 


-1153.70 


5.5007 


Xe 


38.2 


7.6 x 10" 8 


301.1 


Ar 


-369.90 


4.1862 


CH 4 


31.3 


6.0 x 10" 8 


259.1 


Kr 


-603.46 


5.1060 


Kr 


29.2 


5.4 x 10~ 8 


244.4 


Xe 


-819.28 


4.9881 


CO 


25.3 


4.5 x 10" 8 


213.9 








Ar 


22.2 


3.6 x 10" 8 


185.6 








N 2 


21.8 


3.5 x 10" 8 


181.7 
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Figure 1. Stability curves for hydrates, clathrate hydrates (solid 
lines) and pure condensates (dotted lines), together with evolu- 
tionary tracks for the nominal protoplanetary disk model at he- 
liocentric distances of 5, 15 and 20 AU. These evolutionary tracks 
correspond to the thermodynamic pathway followed by the disk 
during its cooling at the considered distances to the Sun (the evo- 
lution of the disk proceeds from high to low temperatures) . The 
abundances of elements are solar, with C0 2 :CO:CH4 = 1:1:1 and 
N 2 :NH3 = 1:1 in the vapor phase. Species remain in the vapor 
phase as long as they stay in the domains located above the curves 
of stability. The label -5.75H 2 designates the clathrate hydrate 
of species X and NH3-H 2 corresponds to ammonia hydrate. 



2.3 Composition of icy planetesimals 

Using the trapping/ formation conditions of the different ices 
calculated at a given heliocentric distance in the outer neb- 
ula (e.g. Table [4]), and knowing their gas- phase abundances, 
one can estimate their mass ratios with respect to H 2 in 
the accreting planetesimals. Indeed, the volatile, i, to wa- 
ter mass ratio in these planetesimals is determined by the 
relation given by Mousis & Gautier (2004): 



Xi 



Xh 2 o Ti(R;Th 2 o, Ph 2 o)' 



(1) 



where Xi and Xh 2 o are the mass mixing ratios of the 
volatile i and H 2 with respect to H 2 in the Solar Neb- 



ula, respectively. £(i?; Ti, Pi) and £(-R; Th 2 o,Ph 2 o) are the 
surface density of the nebula at a distance R from the Sun at 
the epoch of hydratation or clathration of the species i, and 
at the epoch of condensation of water, respectively. From 
rrii, it is possible to determine the mass fraction of species i 
with respect to all the other volatile species taking part to 
the formation of an icy solid: 



Mi 



in, 



j — l,n 



(2) 



with Mi = 1. 



It is important to mention that, after having performed 
calculations using a wide range of values from 10 -3 to 10 _1 
for the viscosity parameter a of our turbulent model, we 
found that the volatile trapping conditions (temperature 
and pressure) remain almost constant at a given heliocen- 
tric distance in the Solar Nebula, whatever the adopted 
value of this parameter. In addition, Alibert et al. (2005a) 
demonstrated that the composition of icy planetesimals 
remains similar whatever their formation distance within 
the same disk, as long as a homogeneous gas-phase is 
postulated within the nebula. These statements imply that, 
whatever the input parameters adopted when modelling the 
disk, and regardless of the formation location considered 
for icy planetesimals at distances beyond 5 AU, their 
composition (in wt%) remains almost constant, provided 
that the gas-phase abundances are homogeneous in the 
nebula. 

Figure [5] shows the variation of the composition of icy 
planetesimals formed in the outer Solar Nebula as a func- 
tion of the C0 2 :CO ratio (shown for values between 1:1 and 
4:1), with CO:CH 4 = 1:1 and N 2 :NH 3 = 1:1 in the gas-phase. 
The composition of ices given in this figure is then valid for 
solids formed at any distance within the outer Solar Nebula, 
provided that the gas-phase is homogeneous and the disk is 
initially warm enough at that location to vaporize the ices 
falling in from the ISM. It can be seen that water remains 
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Figure 2. Plot showing the composition of the ices (wt%) in- 
corporated in planetesimals produced in the Solar Nebula as a 
function of the adopted C02:CO ratio in the initial vapor phase. 
The composition is expressed as a mass ratio (mass of ice i to the 
global mass of ices). The abundances of the considered elements 
are solar with CO:CH4 = 1:1 and N2:NH3 = 1:1 in the vapor 
phase. Note that the mass fractions of Kr and Xe are too low to 
be represented given the scale adopted in the figure. 



Table 5. Ratio of the mass of ice i to the global mass of ices 
(wt%) in planetesimals formed in the outer Solar Nebula, calcu- 
lated for C02:CO = 1:1 and 4:1 in the vapor phase. Both ratios 
are calculated with CO:CH4 = 1:1 and N2:NH3 = 1:1 in the gas- 
phase of the nebula. 



Species 


C0 2 :CO = 


1:1 CC-2 


CO = 4:1 


H 2 


48.8 




25.5 


co 2 


27.5 




58.4 


CO 


7.9 




4.2 


CH 4 


7.6 




3.7 


NH 3 




~ 3.3 




H 2 S 




~ 2.9 




N 2 




~ 1.9 




Ar 




< 0.1 




Kr 




< 0.1 




Xe 




< 0.1 





the most abundant ice in mass, provided that the C02:CO 
gas-phase ratio in the nebula is between 1 and ~1.7. Inter- 
estingly enough, whatever the assumed C02:CO ratio, C0 2 
remains the main carbon species trapped within planetesi- 
mals. We also note that the relative amounts of H2O and of 
other carbon species decrease as the C02:CO ratio increases. 
In contrast, the mass fractions of NH3, N2, H2S, and of the 
noble gases are only weakly influenced by the variation of 
the C02:CO gas-phase ratio. Table [5] gives the mass ratios 
of the ices for the two extreme cases C02:CO = 1:1 and 4:1. 



ciency of this process is particularly dependent on the gas 
density inside the disk, and on the size and density of the 
particles. Cyr et al. (1998) showed that icy particles pro- 
duced in the outer Solar system may have been incorporated 
in bodies formed in the inner regions. Indeed, taking into 
account gas drag, sedimentation and sublimation, these au- 
thors calculated that particles originating from ~ 5 AU may 
migrate to heliocentric distances of the order of 3 AU due 
to gas drag, before subliming as a consequence of the higher 
temperature and pressure domains encountered in the inner 
regions. However, it is important to note that Supulver & 
Lin (2000), Cuzzi et al. (2005) and Ciesla & Cuzzi (2006) 
were unable to reproduce these results. More recently, MA05 
examined the effect of gas drag on the migration of icy par- 
ticles using a time dependent model of the solar nebul4f|. 
They demonstrated that, after the Solar system has expe- 
rienced a few Myr of cooling, it is possible for icy particles 
from the outer Solar system to drift inward to the present 
location of Ceres (~ 2.7 AU). In such an evolved system, 
because the snow-line has already passed through the Main 
Belt, the volatiles contained in the icy particles can avoid 
sublimation during their migration. 

In this paper, we extend the calculations originally per- 
formed by MA05, in order to take into account the possibil- 
ity of volatile capture directly in the form of pure ice, in case 
there is too little water to trap all volatiles as hydrates and 
clathrate hydrates. Moreover, for the sake of consistency, we 
use the disk model that allows the formation of Jupiter and 
Saturn, leading to planets with internal structures consistent 
with observations (see Alibert et al. 2005b). 

The numerical procedure used to calculate the structure 
and evolution of the disk are entirely detailed in Alibert et 
al. (2005c), and in this work we only give the main points. 
More details can be found in the afore-mentioned paper. 
The surface density E in the disk is calculated by solving 
the diffusion equation: 



dE 
dt 



3d_ 

r dr 



1/2 9 , i/ 2 . 
or 



+ E ro (r), 



(3) 



where v is the mean viscosity (integrated along the z axis), 
and Em(r) is the photoevaporation term, taken as in Ve- 
ras & Armitage (2004). The viscosity is calculated in the 
framework of the a-formalism (Shakura & Sunyaev 1973) 
after first solving for the vertical structure of the disk (see 
Alibert et al. 2005c). This calculation also gives us the evo- 
lution of the thermodynamic conditions inside the disk, as a 
function of time and heliocentric distance. These conditions 
are used, together with information on the composition of 
the nebula's gas, in order to derive the amount of volatile 
species trapped inside icy planetesimals. For the model we 
use here, the a parameter is equal to 2 x 10~ 3 , and the total 
evaporation rate is of the order of 10 _8 Mo/year. At the be- 
ginning of the calculation, the gas surface density is given by 
a power law, E oc r~ 3//2 , normalized to have E = 600g/cm 2 
at the current day position of Jupiter. As stated above, this 



3 DELIVERY OF ICY PARTICLES TO THE 
MAIN BELT 

Particles in protoplanetary disks undergo orbital decay due 
to the effects of gas drag (Weidenschilling 1977). The efn- 



5 The work of Cyr et al. (1998) is based on the static Solar Neb- 
ula models elaborated by Cassen (1994). In each of these models, 
the position of the condensation front of water (the "snow-line" ) 
remains at a constant heliocentric distance. In contrast, evolu- 
tionary models cool with time and allow the snow-line to migrate 
closer to the Sun. 
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Figure 3. Thermodynamic conditions in our disk model. Top 
left: midplane temperature profile, top right: midplane pressure 
profile, bottom left: surface density profile, bottom right: aspect 
ratio H/ a profile (H is the semi thickness and a the heliocentric 
distance). The epochs shown in the four panels are, from top to 
bottom, 10 4 yr, 0.5, 1, 1.5 and 2 Myr. 

disk model was used in Alibert et al. (2005b) in order to 
calculate formation models of Jupiter and Saturn. The re- 
sulting thermodynamic conditions are plotted in Fig. [3] for 
early epochs during the evolution of the Solar Nebula. 

Figure U shows the location of the snow- and the 
25K-lintd as a function of time. In our calculations, we as- 
sume that when a particle reaches the snow-line (or the 25K- 
line), it loses all its volatiles (or loses all the volatile species 
whose condensation temperatures are lower than or equal to 
25 K). This is of course a simplification, and more accurate 
calculations should take into account the progressive heat 
diffusion, and the resulting sublimation, within icy parti- 
cles entering the warmer regions of the early Solar Nebula. 
That such low temperatures are reached at small heliocen- 
tric distances in the Solar Nebula can be surprising, given 
that the equilibrium radiation temperatures are of ~160 K 
at 3 AU. However, our calculations are justified by the fact 
that, during the lifetime of the nebula, the Sun's radiation is 
significantly attenuated, even at distances of a few AU, due 
to Rayleigh scattering from molecular hydrogen and dust 
opacity, in particular if no inner gap is postulated in the disk 
(Mousis et al. 2007). In these conditions, light becomes ex- 
tinguished close to the star as a result of the high gas density, 
while the outer regions play little role in the extinction. For 
temperatures below 1500 K, the dominant dimming effect in 
the nebula at wavelengths shorter than a few /im is Rayleigh 

6 The sublimation temperature of water ice is taken as 150K, and 
does not depend on pressure in our calculations. 

7 At temperatures lower than ~25 K, all the main volatile species 
are presumed to be trapped as clathrate hydrates or in the form 
of pure condensates in the Solar Nebula (e.g. Fig. 1). 



scattering from molecular hydrogen (Mayer & Duschl 2005). 
This condition is fulfilled for the entire nebula after 10 5 yr, 
with the temperature beyond 0.6 AU already falling below 
1000 K at this early stage. Moreover, the mass absorption 
coefficient has been estimated to be a m = 5x 10 _4 (cm 2 /g) in 
the Solar Nebula (Mousis et al. 2007). Integrating along the 
mid-plane radial axis of the model from 0.25 AU (the inner 
edge of the disk), the optical depth at 3 AU is initially about 
t = 18. It decreases over time to r = 5.3 after ~5 Myr of the 
disk's evolution. As a result, only ~0.5% of the star's radia- 
tion will be available at this distance within the disk at this 
epoch. We can therefore assume that the disk remains opti- 
cally thick to Solar radiation during at least the first 5 Myr 
of its evolution. Note also that it has been recently showed 
that disks with puffed up inner edge could shield part of their 
outer edge from star light, in agreement with the structure 
and the overall spectral energy distribution (SED) of proto- 
planetary disks around Herbig Ae/Be stars (Dullemond & 
Dominik 2004). This is the case for our solar nebula model 
which is self-shadowed at heliocentric distances higher than 
~2-3 AU (see semi thickness over heliocentric distance pro- 
file in Fig. [3}. Since the Rayleigh scattering, as well as the 
strong opacity of grains impede the light to go through the 
upper layers and illuminate the outer parts of the disk, the 
cooling of the nebula down to low temperatures should fol- 
low, even in the formation zone of asteroids. 

In our calculations, we have considered particles of var- 
ious sizes, ranging from 1 cm to 100 m. We assume that 
they are formed at different locations in the nebula, and at 
different epochs, in a similar way to MA05. Figure [S] shows 
the trajectories of particles of size 1 cm, 10 m, 30 m and 100 
m. For particles ranging between a few centimetres and a 
few metres in size, gas drag is so strong that the trajectories 
would be represented as vertical lines. The trajectories of 
the particles are stopped when they cross the iceline. Fig- 
ure [S] shows that, due to gas drag, particles with the sizes 
considered here can drift to the present day location of the 
Main Belt. The innermost location of particles that have not 
suffered any sublimation of ice is therefore given by the loca- 
tion of the snow-line at a given time. Similarly, the 25K-line 
gives, as a function of time, the innermost location of parti- 
cles that have not lost any of their volatiles. The situation is 
more complex for particles massive enough to significantly 
decouple from the gas, which therefore experience minimal, 
if any, orbital decay (see last panel - 100 m size objects 
- in Fig. [5}. Icy material can thus be transported by such 
medium size particles. If they cross the 25K-line, they would 
be expected to lose some fraction of their volatile content. 
However, particles formed even later, after a few Myr of the 
evolution of the Solar Nebula, but before it becomes opti- 
cally thin, do not encounter temperatures higher than 25K. 
They can therefore preserve all of their volatile species from 
vaporization until their capture by forming asteroids. 

On the other hand, even if icy particles may have drifted 
to heliocentric distances low enough to take part to the for- 
mation of the whole Main Belt without losing a non negligi- 
ble fraction of their volatiles, the current distribution of ices 
that actually exists in asteroids has probably been signifi- 
cantly altered during and/or after their formation. Indeed, 
more heat was generated during the accretion of asteroids 
located at low heliocentric distances, as impacts and col- 
lision velocities were greater. As a result, vaporization has 
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occurred, thus changing the composition of ices incorporated 
in inner belt asteroids. Moreover, the equilibrium radiation 
temperatures reached on the surfaces of Main Belt aster- 
oids after the dissipation of the nebula can be too high for 
the long-term stability of water ice, particularly in the case 
of inner belt asteroids (Jewitt et al. 2007). Indeed, even if 
potential internal heating from the decay of Al 26 and other 
isotopes is not considered, radiation equilibrium may have 
led to the vaporisation of the ice content of the nearer as- 
teroids (semi-major axes of ~2 AU), and melted the ice of 
mid-range asteroids situated at ~3 AU. However, it will 
probably not have affected the ice in the interiors of as- 
teroids located further out. Inner and outer asteroids would 
therefore display no detectable hydratation features either 
because the ice was vaporized and dissipated, or because the 
ice never melted and thus did not react with other miner- 
als that would enable its detection (Cyr et al. 1998). On 
the other hand, mid-range asteroids could be expected to 
have undergone sufficient melting that chemical alteration 
of silicates would occur and be detectable (Cyr et al. 1998; 
Jewitt et al. 2007). Water ice could also be present tem- 
porarily on the surface of large hydrated bodies like Ceres, 
while migrating outwards from sub-surface layers or mantle 
of the asteroids, and before it sublimates into space (Fanale 
& Salvail, 1989). 

These considerations imply that the future detection 
and identification of the volatile phases in asteroids could 
give some constraints on the thermodynamic conditions that 
were present within the Solar Nebula during their accretion, 
as well as during the inwards migration of icy planetesimals. 
For example, assuming a cold accretion for asteroids located 
in the outer belt and that the only detected volatile species 
on the surfaces of some of their members are CO2; H2S, 
NH3 and H2O, one may deduce that the temperature of the 
nebula on the migration pathway of icy planetesimals was 
between ~60 and 80 K (e.g. Table H}. However, we note 
that alteration of the volatile phases in asteroids may occur 
after their formation, as a result of catalytic reactions in 
their interiors. Indeed, this mechanism has been proposed to 
explain the current composition of the plumes released by 
the Saturn's moon Enceladus (Matson et al. 2007) and could 
occur at early epochs following the accretion of asteroids due 
to internal heating caused by the decay of 26 Al. Hence, the 
influence of this chemistry within the interior of asteroids 
may also play a role in the composition of the volatile phases 
that potentially exist in some members of the Main Belt. 



4 IMPLICATIONS FOR THE POROSITIES OF 
ASTEROIDS 

This Section is devoted to a discussion of the determination 
of the taxonomical types of the asteroids, and their densities. 
It appears that some of the results in this area, due to their 
low accuracy, provide poor constraints on the internal struc- 
tures of the asteroids in question, and it is therefore difficult 
to infer their composition from these measurements. This is 
particularly true for the M type asteroids which often show 
a puzzlingly high albedo, and a measured bulk-density much 
lower than the densities of their supposed analogue metallic 
meteorites. This implies that our current knowledge of the 
physical properties of asteroids is not in contradiction with 
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Figure 4. Locations of the 25K- and 150K-lines as a function of 
time for our model of the primitive nebula. 
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Figure 5. Trajectories of solid particles in the primitive nebula, 
for different starting location and starting epochs. The initial lo- 
cation of particles is either 5 AU or 10 AU, and the drift starts 
during the first Myr of the Solar Nebula evolution. The trajecto- 
ries are drawn until particles cross the iceline (heavy solid line), 
from which they are assumed to vaporize. The size of particles is 
given in each panel, and their density is equal to 1 g cm -3 . 



the idea that some of them may have preserved substan- 
tial amounts of the ices accreted during their formation. We 
then discuss the potential influence that the incorporation 
of ices in asteroids may have on their resulting densities and 
porosities. 



4.1 Uncertainties in the taxonomical class and 
densities of asteroids 

The computation of porosity necessitates the knowledge of 
both the bulk-density of an asteroid and the bulk-density 
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of its constituent material. With the exception of the few 
'giant' objects or dwarf planets (Ceres, Pallas and Vesta), 
the bulk-density of asteroids is generally smaller than that 
which would be expected from the densities of supposed me- 
teorite analogues. The relationship between meteorite ana- 
logues and sample densities is, however, not straightforward 
and can, in some cases, be misleading. This is particularly 
important for the objects classified as M, P or E-type aster- 
oids in the Tholen classification scheme. 

Recently, a classification based on spectral morphology 
has been proposed by Bus et al. (2002), where these aster- 
oids are related to class X and sub-classes therein. Aster- 
oids within these classes generally display a linear trend in 
their reflectance curve, and often also have high albedo. In 
the case of radar observations, some M-type asteroids also 
showed particularly strong echo signatures, suggesting they 
are essentially made of metal, just like the M-type mete- 
orites. However, this might not be always the case, and some 
X-type asteroids could also be closer to primordial D- or P- 
type objects, and hence have lower densities. In any case, 
their connection to meteorites is either unknown or unclear 
(Clark et al. 2003; Hardersen et al. 2005) and has been under 
debate during the last decade. M-type asteroids were gen- 
erally thought to be the remnants of the metallic cores of 
differentiated asteroids or planetesimals that would have lost 
their silicate mantle, possibly after a disruptive catastrophic 
collision. There is still debate about the composition and ori- 
gin of M- or X-type asteroids (Busarev 1998; Rivkin et al. 
2000). Following Lupishko & Belskaya (1990), a purely met- 
alic composition is hard to believe, and E-type chondrites or 
stony-iron meteorites could be better analogues. Nowadays 
it appears likely that these classes comprise an assemblage 
of several possibly unrelated surface types with a featureless 
spectrum, not all of which are metallic. Indeed, some aster- 
oids previously classified as M-type have shown spectra with 
the 3 pm signature associated with aqueous alteration which 
clearly makes this metallic-core hypothesis less viable. Sim- 
ilarly, 21 Lutetia has been classified as an M-type asteroid 
since it's albedo is relatively high, but recent spectroscopic 
observations suggest it is actually a classical C-type body 
(Lazzarin et al. 2004). 

On the other hand, the bulk-density of an asteroid 
is derived from measurements of its mass and size. Mass 
determinations obtained as a result of the dynamical per- 
turbations experienced during close encounters are less ac- 
curate, particularly when they are based on a single tar- 
get perturbation (Hilton 2002). One can find large discrep- 
ancies when comparing the values obtained by this tech- 
nique to the one derived from the analysis of a moonlet's 
orbit. For example, the masses of 22 Kalliope and of 87 
Sylvia differ by a factor 2 between Marchis et al. (2003) 
and Kochetova (2004), or between the results of Marchis 
et al. (2005) and Ivantsov (2007), where the latter val- 
ues should theoretically be the less accurate ones. There 
are also large discrepancies between different results ob- 
tained from the same technique but involving different ob- 
servations: depending on the authors, the mass of Psyche 
ranges from (0.68 ± 0.14) x 10" 11 to (0.87 ± 0.26) x lO -11 
or (1.49 ± 0.31) x lO~ n M (Lupishko 2006) or even up to 
(3.38±0.28) x 10 _n or (4.0±1.4) x 10" 11 M Q (Kuzmanoski 
& Kovacevic 2002; Ivantsov 2007), all values that are gen- 
erally outside the usual statistical 3-sigma margin. Finally, 



an error in the size of the body (mainly resulting from the 
determination of the object's albedo) will introduce a large 
error in the determination of the bulk-density. We find that 
the revised density of Psyche from the polarimetric albedo 
of Lupishko (2006), together with the largest value for the 
mass of Kuzmanoski & Kovacevic (2002), would then reach 
the unrealistic value of p = 13.2g/cm 3 . Nonetheless it is 
stressed that the matching between the IRAS-based diam- 
eters and the sizes measured from resolved asteroids in the 
Main Belt is generally good (Cellino et al. 2003; Marchis et 
al. 2006), in contradiction with the latter disagreement of 
Lupishko (2006) for this particular M-type asteroid. 

It therefore seems that it is difficult to obtain mean- 
ingful statistics using the taxonomic M-type class because 
diversity is still present (the situation is barely improved 
when one considers Bus et al.s X-class). It is therefore clear 
that an average value for the density of the M-type asteroids 
(e.g. Krasinsky et al. 2002) can be misleading. It would also 
be difficult to derive knowledge of the internal structure of 
these bodies since the bulk densities quoted in the literature 
can often be in error. 

4.2 Influence of the volatile content in asteroids 
on their bulk densities 

Considering the delivery scenario detailed in section three, 
the asteroids would have been accreted from a mix of icy 
solids that formed initially in the outer Solar Nebula, and 
that preserved their volatile content during their inward mi- 
gration, with heavier and possibly differentiated material 
produced at smaller heliocentric distances. Such a scenario 
would not, however, be applicable to the gravitationally re- 
accumulated aggregates that result from a catastrophic col- 
lision (such as rubble-pile asteroids, and possibly binary ob- 
jects), because their formative collisions would happen at a 
later stage of the evolution of the solar system. Nevertheless, 
ices could still be present in rubble-pile asteroids even after 
exposure to solar irradiation. 

The time scale for re-accumulation scales as (Gp) -1 ' 2 , 
where G is the gravitational constant and p the density 
of the body, so that the dynamical time scale for the re- 
accumulation of the aggregates is - depending on their ac- 
tual velocity dispersion - of the order of several days (Michel 
et al 2001, 2004; Durda et al. 2004; Nesvorny et al. 2006). 
This corresponds roughly to a few revolutions of the debris 
(the remaining material would be ejected). At times larger 
than several weeks, debris are no longer accumulated, but 
rather dispersed into a dynamical family. The time scale 
for re-accumulation after a catastrophic disruption of an as- 
teroid is thus of the order of weeks. On the other hand, 
the sublimation rate experienced by the debris depends on 
several conditions: temperature (heliocentric distance), and 
the presence and thickness of soil and regolith (Chevrier 
et al. 2007). The sublimation rate for a planar surface of 
pure ice is proportional to P s .(M/2ttRT)~ 1 / 2 in a vacuum 
(e.g. Novikov & Vagner 1969, Patashnik & Rupprecht 1975; 
Farmer 1976), where P s is the saturation vapour pressure at 
T, M the molecular weight, R the universal gas constant, 
and T the temperature. 

Depending on the author, the coefficient factor used to 
calculate the sublimation rate from this proportionality in- 
corporates different effects (such as the rugosity, presence of 
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other gas, soil shell burying the ice, etc) and can vary within 
two order of magnitude. The evaporation rate is, however, 
of the order of one metre/year at most (Andreas 2007; Hsieh 
& Jewitt 2006). At 160 K, which is a typical temperature for 
a moderate albedo body in the Main Belt, the recession rate 
of water ice can be as low as a few mm/year (Hsieh et al. 
2004). The sublimation life-time for kilometre sized bodies 
is therefore of the order of a few thousand years. 

The dynamical time for re-accumulation is therefore 
much shorter than that for the sublimation of all the ice in 
collisional debris. As a consequence, one can neglect for our 
purpose the evaporation of ices during the post-impact re- 
accumulation phase. We argue that water and other volatiles 
are not easily revealed, and that such hidden ices will have 
a significant influence on the internal structure and porosity 
of the asteroids. Such a scenario was suggested by Veverka et 
al. (1997) for the asteroid 253 Mathilde, but rejected by the 
authors in the absence of any alteration feature in the spec- 
tra. Similarly, Wilson et al. (1999) consider the case where all 
volatiles must have been completely lost, thus increasing the 
porosity, but again only for the altered bodies. In general, 
the presence of an unrevealed icy fraction will reduce the 
expected bulk-density (the body's average mass-to-volume 
ratio), regardless of its mass distribution, as would do the 
presence of voids and cracks, or any other light material 
such as a regolith. Thus, a low measured density could be 
explained by voids or cracks in a rubble-pile structure, or a 
fractured body, or by the presence of light ices, or both, since 
these hypotheses are not mutually exclusive. In the following 
discussion, we address the bias introduced by the presence 
of ices in a simple computation of the macro-porosity of an 
object. 

The macro-porosity of an asteroid is calculated by the 
determination of the fraction of empty space within the to- 
tal volume, which translates to 1 — pb/p where pb and p 
are the measured bulk-density of the asteroid and the as- 
sumed sample bulk-density, respectively. Following Britt et 
al. (2002), the (uncompressed) sample density is derived 
through knowledge of the asteroid's taxonomic type, using 
values measured in meteoritic analogues, and taking into ac- 
count microporosities in the meteoritic grain. This macro- 
porosity can now be computed by taking into account a pos- 
sible volume fraction of ices ft or silicates f s — V S /(V S + Vi) 
(where V 3 and Vi are the volume of silicates and ices, respec- 
tively), so that f s + fi — 1. This corresponds to an effective 
sample density p s — f 3 (p 3 — pi) + pi, and the bulk-porosity 
is now given by: 



fi (Pi - Ps) + Ps 

where one finds the limiting values of porosity for fi — 
and fi = 1. 

As can be seen in Fig. |6j the bulk porosity of an ob- 
ject with a significant icy component can remain relatively 
low, in the fragmented regime of Britt et al. (2002). The 
corresponding mass fraction is Xi = pi/p s • fi/0- — fi) ~ 
(3 — l))" 1 for typical values of ice and silicates. Thus 
the value of Xi = 17 — 27% for the water content by mass 
determined by McCord & Sotin (2005) in the case of Ceres 
translates to /, « 40 % in volume. Note that the values 
shown in the figure are not strictly restricted to water ice. 
Indeed, the densities of the considered hydrates, clathrate 
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Figure 6. The porosity of asteroids considering a possible 
volume-fraction of ices (/;) or silicates (f s ) of density p = 0.97 
or 3.4, respectively. The ordinate is the average sample density 
(p) = fi {pi — p) + p. The different curves correspond to various 
values of measured bulk-densities (pi,) of asteroids. Upper panel: 
porosity for M-type asteroids with both metal (p = 7.4) and sili- 
caceous components. Bottom panel: porosity for M-type asteroids 
with both metal and ices. 



hydrates or pure ices are very close to that of water ice 
(~0.9 g/cm 3 ). Nevertheless, if a large fraction of ice must 
be present in a given body, it is most likely that they would 
consist of H2O and CO2, as shown in Tabled] 

Assuming a mass fraction of silicate or ices together 
with metallic components would lead to a reduction of the 
mean sample density (p), and hence the computed bulk- 
porosity. As shown in Fig. HJ] the porosity of M-type aster- 
oids 22 Kalliope and 16 Psyche (both X-type in the Bus et 
al. classification) with the adopted updated bulk-densities 
of p = 2.03 g/cm 3 (Marchis et al. 2003) and p = 3.3 g/cm 3 
(Lupishko 2006), respectively, could be reduced by consider- 
ing a substantial volume-fraction of either silicates (olivine, 
pyroxene) or ices (water, ...). In the case of Kalliope, how- 
ever, only a large mass fraction of silicates or ices could 
reduce the porosity to a realistic packing for a supposedly 
post-collisional gravitationally bound granular material (i.e. 
less than approximately 50%). Psyche is a controversial ob- 
ject because its radar echo is high, which makes a metallic 
surface plausible. Keeping in mind that the error bars are 
so large that its density could be in the range [1.7 — 13.2], 
it is clear that the lowest values are feasible only if a sig- 
nificant fraction of lighter material is introduced. Moreover, 
in this latter case, there is no reason to expect a porosity 
larger than 20-30%, i.e. larger than that assumed for frac- 
tured objects. Furthermore, considering that these bodies' 
main constituent is of a lower density (closer to those of sil- 
icates or stony- iron meteorites ~ 3.3 — 4.5), one sees that 
these could still contain a fraction of the order of 30% ice. 

Similarly, one sees that for the low density objects like 
Eugenia, Mathilde and Sylvia and the C-type asteroids in 
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Figure 7. Same as Fig. \E\ with the volume-fraction of ices (fi) 
and silicates (p = 3.4) or the lightest chondrites (p = 2.2). Upper 
panel: porosity for S-type asteroids with some icy fraction (right 
legend). Bottom panel: porosity of C-type asteroids with ices. 



general, introduction of a significant mass fraction of ices 
would change the computed porosity of these bodies. Note 
that a composition of 30% water ice and 70% CI/CM-type 
chondrites in the Cb-type asteroid 253 Mathilde would im- 
ply an average porosity of approximately 30% for both con- 
stituents, a value which seems plausible; at least for the 
porosity of the ice. Pushing this exercise further, one can 
also derive either the density of the other constituents, by 
assuming a constant ice content and a constant porosity, or 
inversely, derive the mass- or volume-fraction of ices at con- 
stant macro-porosity and a given heavy element (chondrite, 
silicate, stony-iron, etc.) density. For instance, the C-type as- 
teroids Mathilde, Eugenia and the X/P-type Sylvia could all 
have the same density within their heavy constituents, and 
a porosity of 40%. This would then require them to contain 
approximately 5, 15 and 45% water-ice, respectively. 

In conclusion, water and other volatiles could have re- 
mained buried inside the Main Belt asteroids since the for- 
mation of the solar System without being easily detected by 
remote spectroscopic observations. The presence of such ices 
will reduce the value of the bulk density of the body that 
would be expected from its taxonomic type, and, moreover 

- given the small amounts of compaction inside asteroids 

- it will additionally decrease the asteroid's macro-porosity 
(Durham et al. 2005). There would be no a priori preference 
of this scenario for any asteroidal taxonomic type (hydrated 
or not), or collisional evolution (given short re-accumulation 
time-scales). On the other hand the size and mass of the as- 
teroid may play an important role, which is not addressed 
here. Neglecting the possible presence of such ices can thus 
increase the error biases on the computation of the porosity, 
leading to it being systematically overestimated. For small 
amounts of ice (fi ~ 0.2, see Fig. [SJ), the error in the poros- 
ity is Sip ~ fi p b (p - pi)/ [p 2 (1 + fi (pi - p)/p)] which yields 
a bias in the calculated porosity of at least 5 to 15 per- 
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Figure 8. The bias introduced in the computation of the porosity 
(ip) by neglecting a possible presence of ices (Eq. 4). Its resulting 
value is given as a function of the actual volume fraction of ices 
(fi) and the sample density (pb) of the body. It is plotted here 
for two cases of pb (solid and dotted grid). 



cent depending on the actual bulk-density of the asteroid 
(pb), its supposed sample density (p), and obviously on the 
actual volatile component (/,). In addition, the precision of 
asteroids bulk densities measurements is increasing with bet- 
ter knowledge of masses and current high resolution angular 
measurements for the size (not to mention space probes). As 
a consequence, the importance of the error or bias discussed 
here should become more and more important. 



5 SUMMARY AND DISCUSSION 

In order to explain the presence of hydratation and cometary 
features in the Main Belt, we have proposed that asteroids 
incorporated during their formation icy particles formed in 
the outer Solar Nebula. We have then calculated the com- 
position of the ices trapped in these planetesimals formed 
beyond a heliocentric distance of 5 AU in the nebula, in a 
manner consistent with the formation of Jupiter, by assum- 
ing that the gas-phase abundances of all elements, in par- 
ticular that of oxygen, are solar. As a result, we have found 
that the ices being formed in the outer Solar Nebula are 
composed of a mix of clathrate hydrates, hydrates formed 
above 50 K, and pure condensates produced at temperatures 
between ~20 K and ~50 K. We have noted that, whatever 
the input parameters adopted in the modelling of the disk, 
or the formation location considered for icy planetesimals 
at heliocentric distances beyond 5 AU, their composition 
remains almost constant, provided that the gas-phase abun- 
dances are homogeneous in the nebula. We have argued in 
this work that gas-drag is responsible for the inward drift of 
icy particles formed in the outer nebula towards the form- 
ing Main Belt. To support this hypothesis, we have showed 
that, at some epochs of the disk's evolution, some particles 
produced in the outer nebula may drift to the current po- 
sition of the Main Belt without encountering temperature 
and pressure conditions high enough to vaporize the ices 
they contain. 
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The current distribution of ices potentially existing in 
asteroids has probably been deeply altered after their for- 
mation. The effect of solar insolation may have vaporized 
the ice within nearer asteroids (semi-major axes of ~2 AU), 
melted the ice of mid-range asteroids situated at ~3 AU, 
but should not have affected the ice in asteroids located 
at greater heliocentric distances. Inner and outer asteroids 
would therefore display no detectable hydratation features, 
either because the ice was vaporized and dissipated, or be- 
cause the ice never melted and thus did not react with the 
surface minerals to a sufficient extent as to allow detection 
(Cyr et al. 1998). In this context, we have proposed that, 
from the detection and identification of initially buried ices 
revealed by recent impacts on the surfaces of asteroids, it 
could be possible to infer the thermodynamic conditions 
that occurred within the Solar Nebula during the accretion 
of these bodies, as well as during the inwards migration of 
the icy planetesimals which they incorporated. However, this 
statement requires that either no parent body processing or 
modification took place during and after the formation of 
asteroids. For example, we have noted that subsequent al- 
teration of the volatile phases in asteroids may occur due to 
catalytic reactions in their interiors. 

We have also investigated the potential influence that 
the incorporation of ices in asteroids may have on their 
porosities and densities. In particular, we have showed that 
the presence of ices can considerably reduce the value of 
the bulk density of the body, and consequently its macro- 
porosity, that would be expected from a given taxonomic 
type. 

That volatiles were delivered to areas within the ice line 
is clearly beyond doubt. In addition to the gas-drag mecha- 
nism described in this work, it is also likely that a significant 
amount of volatile material was dynamically driven inwards 
in the latter stages of planet formation. We still see the tail 
of this dynamical, chaotic volatile movement today - the 
comets we observe passing through the inner solar system 
are the bearers of ices formed far beyond the snow line, and 
held in deep freeze since the early days. During the latter 
stages of planetary migration, the flux of such objects pass- 
ing through the inner solar system, and hence encountering 
the asteroids, was significantly higher. Of particular inter- 
est, when one considers veneers of volatile material near the 
surface of the asteroids, is the Late Heavy Bombardment. In 
the Nice model, (see e.g. Gomes et al, 2005), vast amounts 
of volatile-rich material is flung inwards from the outer so- 
lar system approximately 700 Million years after its birth. 
This event, caused by the resonant destabilisation of the 
outer solar system, would have coincided with a simultane- 
ous stirring of the asteroid belt, leading to an impact flux 
upon the Earth containing approximately even proportions 
of asteroidal and cometary material. It is clear, though, that 
the Earth would not be the only object to encounter volatiles 
injected in this way, and the possibility of a late- veneer of ice 
arriving in the asteroid belt is surely something which must 
be acknowledged in future work. In addition to this aggres- 
sive and chaotic injection of material, there is also a gentler 
mechanism by which volatiles can be driven inwards as a 
result of planetary migration. As planets migrate, material 
can be trapped in the locations of mean-motion resonances 
(MMR), which sweep in front of the planet through it's mo- 
tion. Evidence of material being swept outwards in the reso- 



nances of Neptune is clear for all to see - the Plutino family 
of objects are locked in the 2:3 MMR with the planet, and 
have an inclination distribution which can tell us a great deal 
about the distance over which the planet migrated, sweep- 
ing them along. Inward migration can have the same effect 
- the interior resonances of a planet can collect material as 
it moves inwards, and sweep it along - giving a mechanism 
by which volatile material can be eased inwards, with the 
migration of a giant such as Jupiter. Work such as Fogg & 
Nelson (2006) has shown that such resonant forcing can op- 
erate with a resonable efficiency, even for significantly faster 
migration than expected in our solar system, and so the ef- 
fects of this behaviour should not be ignored in future work. 

In spite of the growing pool of evidence pointing to- 
wards the existence of water ice in the Main Belt, its de- 
tection on asteroids is a challenging observational problem. 
Large bodies such as Ceres are suspected to have retained 
a large amount of water since their formation, perhaps even 
including an internal liquid ocean, throughout the age of the 
solar system (McCord & Sotin 2005). This could particularly 
be the case if this internal water was originally mixed with 
some ammonia, in agreement with our composition calcula- 
tions in section 12.31 which would have the effect of lowering 
the melting point of the water (ammonium bearing minerals 
have been suggested by King et al. (1992) as an alternate 
explanation for the origin of the 2.07 fim band seen in the 
spectrum of Ceres). Nevertheless, internal water can only be 
indirectly probed, either by measuring the hydrostatic shape 
of the object, as was done in the case of Ceres, or by inferring 
its density from its size and mass, when they are known or, 
more evidently, from outgasing activities. The case of Ceres 
is particularly interesting since, in spite of several possible 
pieces of evidence which support it being a highly hydrated 
body, the only report of water detection on the dwarf planet 
is the observation of OH escaping from its northern pol^fl, 
is still not confirmed. Nevertheless, this detection could be 
explained in the context of the accumulation of ice during 
winter on the surface or within the subsurface layer, which 
would then dissipate during summer, when the surface tem- 
perature rises. Similar transient events have been suggested 
as possible mechanisms to trigger the geyser-like activity 
taking place near the south pole of Enceladus and reported 
by Cassini (Porco & Team 2006). 

It is interesting to note that, considering the gravity 
and the day-side temperature of Ceres, any outgassed at- 
mosphere would be rapidly lost. The mean thermal velocity 
vq of H2O, for instance, would be close to the escape veloc- 
ity (voo = 516 km s _1 ). Assuming a subsolar temperature 
of 215 K (Dotto et al. 2000), v would vary between 450 
to 350 km s" 1 from the subsolar point to a zenith angle of 
80° . As a consequence, hydrodynamical escape would occur 
(Vesc/vo ^ 2). The photolysis of H2O by solar EUV makes 
this atmospheric escape even more efficient by giving the 
photodissociation products OH and H some additional ki- 
netic energy. Considering the short lifetime of H2O at ~3 AU 
(<9 days), and the fact that the mean thermal velocity of H 



8 An OH atmosphere was indeed observed around Ceres after 
perihelion by A'Hearn and Feldman (1992) by performing IUE 
long exposure spectra, with column densities of the order of 
10 11 cm~ 2 . 
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atoms exceeds , a tenuous atmosphere of OH is expected 
if water is outgassed by the asteroid at a sufficient rate. Due 
to the transcient nature of the atmosphere, the loss of wa- 
ter to space is limited by the flux of water from the interior 
to the surface. At low latitude, where ice is not stable, the 
continuous flux of water from the interior to space is too 
low to be detected. Only an accumulation of water ice at 
high latitude before perihelion, followed by an outgassing of 
H2O associated with post-perihelion warming seems to re- 
sult in an observable column density of OH. These results 
were found to be consistent with an earlier work done by 
Fanale and Salvail (1989), who estimated the mean loss rate 
of H2O to be in the range 30-300 g s _1 . Even if one assumes 
that the atmospheric loss observed by A'Hearn and Feld- 
man (1992) occurs continuously at the same rate and at all 
latitudes (which is obviously wrong as this maximum loss 
requires high latitudes and post-perihelion conditions), the 
water loss remains below 4 kg s _1 , which, integrated over 
4.5 Gyr, corresponds to only 0.07% of the mass of Ceres. If 
the loss rate of H2O in Ceres remained constant throughout 
its thermal history, the initial water reservoir is thus likely 
to be integraly preserved. Moreover, since the other volatile 
species are expected to be trapped as hydrates, clathrate 
hydrates and pure condensates in this reservoir, we can con- 
clude that they have also been preserved from outgassing 
throughout the thermal history of the asteroid. 

Ceres being the largest and, due to its size, probably 
the wettest Main Belt asteroid, it is an ideal target for 
carrying out observations aiming at constraining its water 
regime. The experiment searching for water being vaporised 
within the polar regions of Ceres should be repeated with 
the state-of-the-art instrumentation available today on large 
telescopes. Such a detection would confirm unambiguously 
the presence of a large amount of water near the surface of 
Ceres. Direct observation of water ice, or of the effects of 
hydration, on the surface of Ceres can also be attempted 
for lower latitudes on the asteroid using a combination of 
high-angular resolution and spectroscopic instruments per- 
mitting the full resolution of its surface to the 30-40 km 
level. Due to its low spectral resolution, imaging of the sur- 
face of Ceres, even when it is spatially resolved using HST 
or adaptive optics, is not sensitive to the presence of ice, 
while the detection of such is within the reach of low resolu- 
tion spectroscopic observations (e.g. the detection of absorp- 
tion features in the 1.0-3.5 /xm region). A spatially resolved 
spectroscopic mapping of the surface of Ceres in the near- 
infrared can be done with today's ground-based telescopes 
and would permit the mapping of the strength of the 3 fj,m 
band, and allow the search for regions on the surface where 
interstitial water ice, or hydration features could be present, 
for instance at the location of cracks within the surface of 
Ceres, or the locations of deep impact craters. Indeed, recent 
HST (Thomas et al. 2005) and adaptive optics (Carry et al. 
2007) imaging observations of Ceres revealed the presence 
of large impact craters across its surface which have likely 
disrupted the outer crust of the asteroid enough to directly 
expose the sub-surface mantle of wetter material. Finally, 
a spectroscopic study of the surface of Ceres, in order to 
search for the spectral signature of water and maybe those 
of other volatiles, should not be limited to one wavelength 
region (although the near-infrared range offers many diag- 
nostic bands) but should, instead, encompass a wider range, 



from the near-UV to infrared wavelengths, in order to im- 
prove the identification of the chemicals species responsible 
for these spectral features. 

Finally, the NASA Discovery mission Dawn, which has 
been launched in September 2007 and whose arrival at Ceres 
is scheduled for 2015, will certainly bring new constraints on 
the presence of volatiles in the Main Belt. In particular, the 
Dawn mapping spectrometer (MS) covers the spectral range 
from the near UV (0.25 u.m) through the near IR (5 u.m) 
and has moderate to high spectral resolution and imaging 
capabilities (Russell et al. 2004) . These characteristics make 
it an appropriate instrument for determining the asteroid's 
global surface composition. Near-infrared mapping of the 
surface of Ceres at small spatial scales will be very sensi- 
tive to volatile concentrations and may reveal ice spots on 
fresh impact-crater ridges. Moreover, the gravitation inves- 
tigation of Ceres will allow the determination of its gravity 
field up to the 12th harmonic degree (Russell et al. 2004). 
Such a measurement will enable the shape and gravity mod- 
els to characterize crustal and mantle density variations and, 
consequently, the amount of volatiles trapped therein. 
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